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Abstract 

Using a density functional based interface displacement model we determine 
the effective interaction potential between two spherical particles which are 
immersed in a homogeneous fluid such as the vapor phase of a one-component 
substance or the A-rich liquid phase of a binary liquid mixture composed of 
A and B particles. If this solvent is thermodynamically close to a first-order 
fluid-fluid phase transition, the spheres are covered with wetting films of the 
incipient bulk phase, i.e., the liquid phase or the B-rich liquid, respectively. 
Below a critical distance between the spheres their wetting films snap to a 
bridgelike configuration. We determine phase diagrams for this morphological 
transition and analyze its repercussions on the effective interaction potential. 
Our results are accessible to force microscopy and may be relevant to floccu- 
lation in colloidal suspensions. 
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I. INTRODUCTION 



In view of understanding a particular phenomenon in condensed matter, theory is sup- 
posed to identify the corresponding relevant degrees of freedom and to provide the effective 
interaction between them by, approximately, integrating out the remaining ones so that one 
is left with a manageable model. It is a major challenge to determine the effective interac- 
tions because that requires to calculate the partition function of the whole system under the 
constraint of a fixed configuration of the relevant degrees of freedom. The benefit for carry- 
ing out this constrained calculation, which in general is more difficult than the original full 
problem, is twofold. First, there is a gain in transparency by describing the system in terms 
of relevant degrees of freedom. Secondly, it is typically less risky to apply approximations 
for the partial trace because they only concern the less relevant degrees of freedom. 

The determination of the phase behavior and of the structural properties of multi- 
component fluids represents a case study for this general approach. If the composing par- 
ticles of the mixture are of comparable size and shape their degrees of freedom have to be 
treated on equal footing. The well developed machinery of liquid state theory JE] offers var- 
ious techniques to cope with this problem. However, these techniques fail to yield reliable 
results if, e.g., one component is much larger than the others; in this case numerical sim- 
ulations become inefficient and integral theories lose their accuracy. Colloidal suspensions 
are a paradigmatic case for such highly asymmetric solutions. For their description these 
difficulties can be overcome by resorting to the general scheme laid out at the beginning with 
the positions of the colloidal particles as the relevant degrees of freedom. Accordingly the 
degrees of freedom of the small solvent particles are to be integrated out for a fixed config- 
uration of the colloidal particles which we assume to be smooth, monodisperse spheres. At 
sufficiently low concentrations of the suspended particles this leads to an effective pair po- 
tential between them. In many cases the effective potential resembles the bare one, i.e., the 
one in the absence of the solvent, but with modified, effective interaction parameters which 
depend on the thermodynamic variables of the system such as pressure and temperature. 

The effective pair potential acquires additional new features if the solvent is enriched with 
particles of medium size such as, e.g., polymers. If the colloidal particles come close to each 
other the depletion zones around them, generated by the finite size of the medium particles, 
overlap leading to an entropically driven attraction of the colloidal particles [|2],|3| . Correlation 
effects can modify the form and the range of these depletion forces considerably ||,[5|]. These 
effective potentials have indeed turned out to be successful in describing the phase behavior 
of colloidal suspensions 0. 

Qualitatively new aspects arise if the solvent particles exhibit a strong cooperative be- 
havior of their own such as a phase transition which proliferates to the effective potential 
between the large particles. If the solvent undergoes a continuous phase transition, thermal 
Casimir forces between the large particles are induced due to the geometrical constraint 
they pose for the critical fluctuations f7],f|. Such forces are long-ranged and have a strong 
influence on the phase behavior of the colloidal particles P,nj. If the solvent is thermody- 



namically close to a first-order phase transition, wetting phenomena JTT]] can occur at the 
surfaces of the dissolved particles (see Ref. [|l^] and references therein). If the bulk phase 
of the solvent is the vapor phase of a one-component fluid, the surfaces of the large spheres 
can be covered by a liquidlike wetting film. This situation corresponds to aerosol particles 
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floating in a vapor. If the bulk phase of the solvent is the A-rich liquid phase of a binary 
liquid mixture composed of (small) A and B molecules, the dissolved colloidal particles can 
be coated by the B-rich liquid phase of the mixture. If the wet spheres approach each other, 
at a critical distance the two wetting films snap to a bridgelike structure. This morphological 
phase transition is expected to yield a nonanalytic form of the effective interaction potential 
between the large spheres. This nonanalyticity demonstrates that cooperative phenomena 
among those degrees of freedom which are integrated out can leave clearly visible finger- 
prints on the effective interaction between the remaining relevant degrees of freedom. The 
study of this kind of profileration is not only of theoretical interest in its own right but 
seems to play an important (albeit not exclusive |13[) role for the experimentally observed 
flocculation of colloidal particles dissolved in a binary liquid mixture close to its demixing 
transition into an A-rich and a B-rich liquid phase P^I5|. This observation has triggered 
numerous theoretical efforts devoted to various possible explanations of it. Since they are 
reviewed in Sec. I of Ref. |T2j and more recently in Ref. |TP[ the interested reader is referred 
to there and we refrain from repeating this discussion here. 

In our present analysis of this problem we apply density functional theory |2(J which offers 
two advantages. First, this technique is particularly well suited to calculate, as required 
here, free energies under constraints. Secondly, it allows one to keep track of the basic 
molecular interaction potentials of the system. We focus our interest on thermodynamic 
states of the solvent which are sufficiently far away from its critical point so that the emerging 
liquid-vapor interfaces of the wetting films exhibit only a small width. Therefore we can 
apply the so-called sharp-kink approximation which considers only steplike variations of the 
solvent density distribution and thus leaves the interface position as the main statistical 
variable. This approximation has turned out to be surprisingly accurate for the description 
of wetting phenomena [|21||. Our analysis extends and goes beyond previous efforts 



which are based on a similar interface displacement model grounded on a phenomenological 
ansatz. Whereas Refs. and [23] are aimed at mapping out the phase diagram in terms of 
interaction parameters for the bridging transition mentioned above, we focus on the effective 



interaction potentials between the wet spheres, which are not presented in Refs. |22j and |23 
and on their microscopic origin. Inter alia, this allows us to compare the effective interaction 
potential between the colloidal particles with the bare one, i.e., in the absence of the solvent, 
and thus to comment on the quantitative relevance of the solvent-mediated interaction. 
Moreover, we present the phase diagram of the system in terms of the thermodynamic 
variables temperature and chemical potential which is also not contained in Refs. |22[ and 

II- 



In Sec. [II] we describe the implementation of a simple version of density functional theory 
for the present problem. For reasons of simplicity we confine our analysis to liquid-vapor 
coexistence of a one-component solvent; the generalization to a binary solvent is straight- 
forward. In Sec. [TTI| we present some examples for the numerically calculated wetting film 
morphologies and discuss a phase diagram for the aforementioned morphological transition, 
and in Sec. [TV] we analyze the effective wetting-induced interaction potential between the 
spheres as a function of the distance between the spheres and the undersaturation. The ex- 
perimental relevance of our model calculations is discussed in Sec. |V|and Sec. |VT| summarizes 
our main results. The Appendix contains some technical details. 
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II. DENSITY FUNCTIONAL THEORY 



A. Model 



We consider two identical, homogeneous, and smooth spherical particles of radius R 
whose centers of mass are separated by a distance D (see Fig. [3]). They are immersed in a 
fluid of particles of number density p(r) which interact via a Lennard- Jones potential 
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The system is symmetric with respect to a rotation around the axis which connects the 
centers of mass of the spheres (Fig. |I|) and with respect to a reflection at a plane in the 
middle between the spheres that is perpendicular to the symmetry axis. Since we work in a 
grand canonical ensemble and the fluid particles are subject to the external potential exerted 
by the spheres, the equilibrium number density profile of the fluid particles exhibits these 
symmetries, too. Therefore we describe the system in cylindrical coordinates (r±, <p, z), with 
the z axis being the symmetry axis of the system. The two centers of mass of the spheres are 
located at (r = 0, z = ±D/2) such that the spheres occupy the volumes S± = {r(r±, 0, z) = 
(x,y,z) = (V ± cos0,r ± sin0,2) E R 3 \±D/2-R <z< ±D/2 + R, yfi* + {z =f D/2) 2 < R}. 
The external potential exerted by both spheres on each individual fluid particle is 

v ta t(r±, z-R) = v (^r\ + (z - D/2) 2 ; R J + v (^Jr 2 ± + {z + D/2) 2 ; R^j (2.2) 

where (see Eq. (A.4) in Ref. pf) 

( m 9 ( 1 1 \ ( 1 1 



8 y V^(^ + ^) 8 r(r-R) 8 J M \(r + R) 9 (r - R) 



2 \r(r + R) 2 r(r - R) 2 J \ (r + ,R) 3 (r - R) 

is the interaction potential between a single sphere of radius R and a fluid particle at a 
distance r > R from the center of mass of the sphere. In a continuum description, v(r; R) 
follows from an integration of the Lennard- Jones potential 

a sf\ 12 (<T,f 



Mr) =4^^-) -^j ) (2.4) 

between a molecule of the spherical substrate and a /luid particle. The subscript sf denotes 
the parameters of the dispersion interaction between a particle in the fluid and a particle 
in the spheres. One has 113 = ^es/p^crfj and u 9 = ^^sfPs^lj where p s is the number 
density of the particles forming the spheres. (Many colloidal particles exhibit an even more 
complicated substrate potential because they are coated by a material different from their 
core so that they are no longer radially homogeneous as assumed for Eq. (|2.3|) .) 

Within our density functional approach the equilibrium particle number density distri- 
bution of the inhomogeneous fluid surrounding the spheres in a grand canonical ensemble 



minimizes the functional 20 
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fi([p(r)];T,p) = J d 3 r(f HS (p(r),T) + (v tot (r)-p)p(r)) 

+\J J d 3 rd 3 r'w{\r-r'\)p(r)p(r'). (2.5) 



/ 



Vf — V \ (S + U5_) is the volume accessible for the fluid particles, V is the total volume 
of the system; V — > M. 3 in the thermodynamic limit. Equation ( |2.5| ) does not include the 
bare interaction potential $(-D; R) (see, c.f., Sec. |V[) between the solid spheres, separated by 
vacuum, generated by the dispersion forces between the molecules forming the two spheres. 
fns(p, T) is the free energy density of a hard-sphere fluid of number density p at temperature 
T. In Eq. (|2.5|) , the hard-sphere reference fluid is treated in local density approximation. 
We apply the Weeks-Chandler-Andersen procedure |24[] to split up 0(r) into an attractive 



part 4>att( r ) and a repulsive part 4> rep (r). The latter gives rise to an effective, temperature 
dependent hard-sphere diameter 

2 1 /6 CT 

d{T)= J dr (l-exp(-^p)) ( 2 -6) 
o 

which is inserted into the Carnahan-Starling expression |25| 

fns( P , T) = k B Tp (ln(pA 3 ) - 1 + j£^f) ^ 

for the free energy density fns of the /iard-sphere fluid, where rj = ^p(d(T)) 3 is the dimen- 
sionless packing fraction and A is the thermal de Broglie wavelength. We approximate the 
attractive part of the interaction 4> a tt(j) by 

w{r ) = ^l {r 2 + a 2 )- 3 (2.8) 



f 32 

^ rw ( r )= d 3 r(f) att {r) = - — v^vrea 3 (2.9) 

JR3 9 



with 

w = 

in order to simplify subsequent analytical calculations. The double integral in Eq. ( |2.5| ) 
takes into account this attractive interaction within mean-field approximation. 

In the bulk the particle density p 7 (where 7 = l,g denotes the liquid and vapor phase, 
respectively) is spatially constant, leading to (see Eq. ( |2.5|) ) 

fi fe (p 7 , T > P) = fHs(p 7 , T) + ^w p 2 - p,p 7 (2.10) 

for the grand canonical free energy density of the 6ulk fluid. Minimization of Qf, with respect 
to p 7 yields the equilibrium densities. The line p = po(T) of bulk liquid- vapor coexistence 
and the two bulk densities pi and p g at coexistence follow from 



dtt h 



dp 



dtt h 



p=p s d P 



and n b {p g ) = Q b { Pl ). (2.11) 
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p=pi 



For p 7^ p , i.e., off coexistence, only the liquid or the vapor phase is stable. In this case the 
density of the metastable phase corresponds to the second local minimum of Q b . 
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B. General expressions for the contributions to the effective interaction potential 



Henceforth we consider the case that the substrate is sufficiently attractive so that the 
liquid phase is preferentially adsorbed. Therefore, if in the bulk the vapor phase is stable 
(/i < /i ), the fluid density is significantly increased in the vicinity of both spheres. In the 
spirit of the so-called sharp-kink approximation (see Sec. | and Ref. pl| ) we assume that a 



thin film of constant density p\ but with locally varying thickness is adsorbed at the surfaces 
of the spheres, separating the spheres from the bulk vapor phase of density p g . This wetting 
film encapsulating both spheres is characterized by a function h(z): 

p(r) = p(r ± , 0, z) = G{r ± -{R + d s ))(e{h{z) - r ± ) Pl + Q{r x - h{z))p g ^ (2.12) 

where denotes the Heaviside step function. The length d s takes into account the excluded 
volume at the surfaces of the spheres which the centers of the fluid particles cannot pene- 
trate due to repulsive forces. The profile h(z) as given by Eq. ( |2.12|) can describe both a 



configuration in which the wetting films surrounding each sphere are connected by a liquid 
bridge as well as the configuration in which both single spheres are surrounded by disjunct 
wetting layers. In the latter configuration there is a region around z = with h(z) =0. 

Inserting p(r±, 0, z) from Eq. (|2.12j ) into the functional Q in Eq. fl2.5|) leads to a decom- 
position of Q = Vol(V/)fi;,(p(,) + Qs into a bulk and subdominant contributions. The bulk 
contribution is Vol(V/)fi;,(p 9 ) (with Qb given by Eq. ( |2.10|) ) and corresponds to the vapor 
phase which is stable in the bulk. The subdominant contribution is 

n s [h] = n s i + n ex [h] + n ei [h] + n lg [h] (2.13) 

where only Q s i is independent of h(z) and all the other three contributions are functionals 
of h(z). Since we have not found an indication for spontaneous symmetry breaking, in the 
following we discuss only symmetric configurations with h(z) = h(—z). 

n ex [h(z)} = voi(£)(p b ( Pl ) - n b ( Pg )) (2.14) 

with 

L z 

Vol(£) = 2ir J dzh 2 {z) - yi? 3 (2.15) 


is an excess contribution which takes into account that the volume C = JC \ (<S_ U«S+) is filled 
with the metastable liquid instead of the vapor phase; K, = {r(r±,(p, z) G R 3 |r^ < h(z)} is 
the volume enclosed by the liquid-vapor interface. (The excluded volume due to d s enters 
into fl s i (see, c.f., Eq. ( [2.21|) ).) This free energy contribution Q ex vanishes at two-phase 
coexistence p = po(T) (compare Eq. ( |2.11| )). 2L Z is the extension of the total volume of the 
system V in 2 direction; L z — > 00 in the thermodynamic limit and h{z > z max ) = with 

Q ei [h{z)} =2Ap J d 3 r( Pl (t(r,S^+t(r,S + ))-v tot (v)^ (2.16) 
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can be interpreted as the integrated effective interaction between the spheres and the liquid- 
vapor interface described by h(z); Ap — pi — p g . V_ is that part of the volume V with z < 
(we note again that V_ — > M?_ in the thermodynamic limit which is always considered here), 
analogously /C_ is the part of the set K with z < 0. In Eq. (|2.16| ) we have introduced the 
interaction potential 



t(r;M) = J dVtoflr-r'l) 
M 



(2.17) 



between a fluid particle at r and a region M. (with r ^ M.) homogeneously filled with the 
same fluid particles (analogous to the function t(z) introduced in Refs. |IT| and |21| in the 
case of a planar substrate). v tot is the total interaction potential between a fluid particle 
and both spheres (see Eq. (|2.2|) ). Finally, 



(2.18) 



Q lg [h(z)] = -(Ap) 2 J d 3 r (*(r;/C_)+*(r;£ + )) 



is the free energy contribution from the free /iquid-gas interface. It is a nonlocal functional 
of h(z) in contrast to Q ex and Vt ei whose dependence on h(z) enters only via the integration 
volume /C_. The local approximation thereof, which is provided by the gradient expansion 
of Eq. gngp , is 



[2.19) 



In Eq. (|2~T9| ) 



a 



oo oo 



(p) 

ig 



-^(Ap) 2 J dz J dz' J d 2 r\\w (yjr\ + z' 2 



^2.20) 



is the interfacial tension of a planar, free liquid- vapor interface in sharp-kink approximation. 
We note that, strictly speaking, the surface tension of a curved liquid- vapor interface depends 
on the local radius of curvature (see Fig. 2 in Ref. and the references therein concerning 
the Tolman length). This curvature dependence is omitted in the local model presented here. 
However, for spheres of radius R > 20a as considered henceforth the curvature correction 
is less than 1%. Similar arguments hold for the deviation of the actual liquidlike density in 
the wetting film from the bulk value pi. 

For our choice of interaction potentials (Eq. ( |2.1| ) and ( |2.4j )) a tedious calculation leads 
to explicit expressions for the contributions Q e i and Qi g which are given in the Appendix. 
The remaining contribution 

Vsi = -Pi / d 3 r (p,(t(r, <S_) + t(r, 5+)) - 2^(r)J - Q b ( Pl ) -^-((R + d s ) 3 - R 3 ), (2.21) 

which is independent of h(z), is the sphere- /iquid interfacial free energy corresponding to the 
interface between the spheres and the liquid phase. The last term in Eq. (|2.21| ) takes into 
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account the excluded volumes at the surfaces of the spheres. In the limit of large separations 
D one has 



Q sl {D -> oo) - 2^ ~ IT 6 (2.22) 

with the sphere-liquid interfacial free energy iV^ of a single sphere immersed in the liq- 
uid phase. The leading power law ~ D~ 6 in Eq. ( |2.22| ) can be inferred from the following 
consideration: if present, the second sphere displaces a spherical volume from the homoge- 
neous liquid phase so that the free energy of the interaction of the first sphere with the bulk 
liquid is reduced by the interaction free energy of that sphere with the displaced spherical 
liquid volume. This latter interaction decays as D~ 6 for large separations D, at which the 
dispersion interaction between two spherical objects resembles the dispersion interaction be- 
tween two pointlike particles. (Here, as before, we have not yet taken into account the bare 
interaction potential $(-D; R) between the two solid spheres; but see, c.f., Sec. |Vj.) 

Up to the bulk contribution the grand canonical potential of the system is the minimum 
of Q s [h(z)} with respect to the profile h(z): 

tt s = tt s (D; R) = min (£l[h(z)]). (2.23) 

{h(z)} 

Thus the equilibrium interface morphology h(z) minimizes fig[/i(z)] which includes the con- 
tributions fl ex [h(z)}, ft ei [h(z)], fli g [h(z)}, and Q s i. The functional used in Refs. [^] and |2B 



(Eq. (1) in both references) is, albeit formulated in another coordinate system and using a 
more phenomenological ansatz for the basic interaction potentials, essentially identical with 
the sum (fl\° c + Q ex + Q e i)[h(z)]. However, this model description does neither incorporate 
the bare dispersion interaction of the two spheres (of., Sec. |V|) nor the free energy contri- 
bution Q s i which describes the sphere-liquid interfacial free energy. We emphasize that the 
consideration of the contribution Q s i - which does not depend on h(z) - is not essential for 
the determination of the equilibrium wetting film morphology and hence it is not relevant for 
the thermodynamic phase diagram of thin-thick and bridging transitions (Fig. 2 in Ref. p3| ) 
for a fixed separation D between the spheres. But the term fl s i is crucial to the shape of the 
effective, wetting-induced interaction potential between the spheres, i.e., its dependence on 
D (see Eq. (ggp ). 



III. MORPHOLOGY OF THE WETTING LAYERS 
A. Interface profiles 

The actual wetting layer morphology h(z) follows from numerical minimization of the 
functional f2s[/i(,2)] (Eq. (|2.13|) ) for a given temperature T and undersaturation A/i = 
Ho(T) -//, with the contributions VL ex (Eq. (|2TJ)), Vt ei (Eq. (pT6|)), Vt sl (Eq. fl£2iP), and 



Qi g (Eq. ( p,18| ) within the nonlocal and Eq. Q2.19| ) for the local theory). Within a range of 
parameters (T, A/i) the numerical minimization yields two different solutions for h(z), one 
with a liquid bridge and one without, depending on the initial function h(z) used in the it- 
eration scheme for the minimization. For small separations a 2R only the solution which 
exhibits a liquid bridge is stable whereas for large separations a ^> 2R only the solution 



S 



without bridge minimizes Qs- For large distances D the minimization consistently yields 
twice the result known for a single individual sphere enclosed by a wetting film (compare 
Ref. ||12|| ). This observation amounts to a useful check of the numerical procedure. 

As a first example, in Fig. |2| we present the numerical results for a wetting layer enclos- 
ing two spheres of radius R = 20a. For our particular choice of interaction parameters, at 
coexistence A/i = the wetting film on each of the single spheres alone exhibits a first- 
order thin- thick transition (which is the remnant of the first-order wetting transition on the 



corresponding planar substrate, see Fig. 8(a) in Ref. |12[) at T* t = kgT/e pa 1.271 (which 
corresponds to T a /T c pa 0.9 where T c is the critical temperature of gas-liquid coexistence in 
the bulk). The planar substrate, i.e., a single sphere in the limit R — > oo, exhibits a genuine 
first-order wetting transition (with the film thickness jumping to a macroscopic value) at 
T* pa 1.053 (T w /T c w 0.75, T tt /T w w 1.21). Figure |(a) depicts a typical solution with a 
bridge, here for a separation a = D — 2R = 10a (D = 50a) and the thermodynamic param- 
eters T* = 1.3 > T£ and A/i = 0, i.e., at liquid- vapor coexistence. The solution without a 
bridge for the same choice of parameters is shown in Fig. 0(b). The latter solution has a 
higher free energy than the former one. Therefore the solution with bridge is thermodynam- 
ically stable whereas the solution without bridge is metastable. For the solution without 
a bridge the distortion of the liquidlike layer around one sphere due to the presence of the 
other sphere is not visible. Finally, Fig. 0(c) displays the wetting film morphology for the 
stable state with bridge at the temperature T* = 1.2, i.e., below the thin-thick transition 
temperature T^. (We note that the thin-thick transition temperature T u for each sphere is 
slightly shifted by the presence of the second sphere. However, as already pointed out in 
Ref. ||23|| , this effect is negligibly small.) In any case, the difference between the nonlocal and 
the local theory is very small. This latter result is in accordance with the findings for the 
comparison between the nonlocal and the local description of the three-phase contact line 
on a homogeneous substrate and of the wetting layer morphology on a chemically structured 
substrate (compare Ref. PB|). For this reason, henceforth we only consider the local theory. 

Figure |3] shows another pertinent example. Here we study the wetting layer morphology 
for two larger spheres of radius R = 50a as a function of the undersaturation A/i along the 
isotherm T* = 1.2. The interaction potential parameters are the same as for the previous 
first example and the separation of the surfaces a is 20a (D = 120a). At coexistence each 
single sphere exhibits a first-order thin-thick transition at T£ ~ 1.191 (i.e., T tt /T c pa 0.84 and 
Ttt/T w pa 1.13). In analogy to the prewetting line on a homogeneous substrate there is a line 
of thin-thick transitions (T, A/i tt (T)) which intersects the liquid-vapor coexistence line at 
(T = T u , A/i = 0) (compare with, c.f., Fig. [|and Fig- 8(a) in Ref. [0). At the temperature 
T* = 1.2 > T* t considered here the thin-thick transition occurs at A/i* t = Afi tt /e pa 0.0103. 
Upon reducing the undersaturation along the isotherm, starting at, e.g., A/i* = 0.05, first 
the configuration with thin films and without bridge is stable (Fig. 0(a)). For A/i < A/i& 
(bridging transition) with A/i* 4 pa 0.0235 > A/i* t (T) the solution with bridge becomes stable, 
but the layers enclosing the spheres still remain thin (Fig. §(b)). Upon further reduction of 
A/i, at Afi tt (T) the second transition from a solution with bridge and thin films to a solution 
with bridge and thick films (Fig. 0(c)) takes place. (As before, concerning the value of T£ t at 
coexistence, also the value A/i* t (T) is practically unchanged by the presence of the second 
sphere - even for the bridge configuration.) We note that for this choice of parameters and 
in the case of a solution with bridge and thin films (Fig. 0(b)) the profile h(z) exhibits six 
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turning points instead of only two as for the case of a solution with bridge and thick films 
(Fig. y(c)). This rich curvature behavior is caused by the details of the effective interaction 
potential between the spherical substrate surfaces and the liquid-vapor interface (see Sec. 2.3 
in Ref. ||12||), similar to the curvature behavior of the liquid- vapor interface when it meets 
a homogeneous, planar substrate forming a three-phase contact line (compare Ref. [f26|j). 
These features may also occur for a bridge configuration with thin films at coexistence and 
T < T tt . 



B. Phase diagram 

The example presented in the previous paragraph shows that besides the gas-liquid co- 
existence curve A/i = the T-A/i phase diagram of the system contains two distinct lines 
of first-order phase transitions: a line of thin-thick transitions (T, Afi tt (T)) on the single 
spheres (which is the remnant of the line of prewetting transitions on the corresponding 
flat substrate and which is, as stated above, practically unshifted by the presence of the 
second sphere) and a second, independent line of bridging transitions (T, A/i^(T)). If one 
crosses the latter along an isotherm T = T approaching coexistence (T ,A/i — > 0), at 
A/i = Afi bt (T ) a transition from the configuration without bridge (A/i > A/i M (T )) to 
a configuration with bridge (A/i < A/i W (T )) occurs. The derivative dfls/dAfi is discon- 
tinuous at A/iM, indicating that the bridging transition is first order. Figure |] shows the 
T-A/x phase diagram for the two spheres with R = 20a for D = 50a (a = 10a). The line 
of thin-thick transitions intersects the liquid- vapor coexistence line at T t * ~ 1.271 with a 
finite, negative slope (compare Fig. 8(a) in Ref. |r2"||). It extends into the vapor phase region 
(A/i > 0) of the phase diagram and ends at a critical point. The line of bridging transitions 
intersects the coexistence line also with a finite, negative slope. On the other end, within 
our sharp-kink interface model, it happens to be cut off at that metastability line in the 
phase diagram at which the second minimum of the bulk free energy at high fluid density 
(Eq. ( [2.10|) ) ceases to exist so that for larger undersaturations the liquid phase is not even 
metastable. Within a more sophisticated approach, e.g., by seeking the full minimal density 
distributions of Eq. ( |2.5| ), the line of bridging transitions is expected to end in a critical 
point, too. (Concerning the effect of fluctuations on these mean field predictions see the 
following paragraph.) The line of bridging transitions is entirely located in the region where 
the liquidlike films on the spheres are thin. Moreover, the effect of the presence of the liquid 
bridge on the line of thin-thick transitions is negligibly small. In Fig. ^ the relative location 
of the bridging transition line and the thin-thick transition line corresponds to our specific 
choice of the interaction potential parameters as well as the chosen size of and distance be- 
tween the spheres. Changing these parameters will lead to shifts of these lines and, possibly, 
to different topologies of the phase diagram. Here we refrain from exhaustingly presenting 
all possibilities which can occur according to Refs. |22| and [ p3| . 

Since the liquid volume enclosed by the interface h(z) is quasi-zerodimensional, fluc- 



tuation effects destroy the sharp first-order phase transition (see Refs. [57| and |28[). In 
Sec. 4 of Ref. [0 it has been extensively discussed how finite size effects smear out the 
thin-thick transition such that the thickness increases sharply but continuously within a 
range <5/i around A/i t4 (T); these results apply analogously to the present case. Using similar 
approximations we obtain a range 5/i between <5/i* ~ 0.004 for T* = 1.16 and 5fi* ~ 0.02 
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for T* = 1.26 over which the bridging transitions shown in Fig. f| are smeared out around 
Afib t (T). Thus close to A/i = the quasi- first-order thin-thick transitions are clearly visi- 
ble. However, for larger values of A/i they become progressively smeared out such that their 
critical points predicted by mean field theory are erased by fluctuations. 



IV. EFFECTIVE FILM-INDUCED INTERACTION POTENTIAL 

A. Shape of the effective potential, metastability, and asymptotic behavior 

In the following we change our point of view: we vary the distance D between the centers 
of mass of the spheres instead of the thermodynamic parameters T and A/i. Figure |5| shows 
the grand canonical potential Qs corresponding to the wetting layer morphologies for the 
case R = 20a and T* = 1.2 (Fig. 0(c)) as a function of the separation a = D — 2R for several 
values of A/i. Qs is the minimum of Qs[h{z)} (Eq. ( |2.13 )) for the given set of parameters 



T, A/i, and D = 2R + a. For each value of A/i there are two branches of the free energy, 
one corresponding to the solution without bridge, which for the case R = 20a considered 
here exists only for a > 0.15.R, and the other corresponding to the solution with bridge 
which exists up to a « 0.65.R and a ~ 0.6R for A/i* = and A/i* = 0.01, respectively. At 
a certain value D = Db t or, equivalently, a — a^t, which are functions of A/i, a first-order 
phase transition occurs with discontinuous derivative dVLs/dD between the solutions with 
and without bridge. The main effect of increasing the undersaturation A/i is that the free- 
energy curves are rigidly shifted upwards. This shift is approximately proportional to A/i 
and larger in the case of the solution with bridge, resulting in the dependence of Du on A/i. 
The values of Qs shown in Fig. are obtained within the local theory. The nonlocal theory 
yields the same functional dependence £ls(D) but with a slight and rigid shift of the free- 
energy curves, relative to the results of the local theory, of the order of 0.1% and of the same 
sign and size for both the solutions with and without bridge. Finite-size effects again destroy 
the sharp first-order bridging transition; we obtain a range 5D ~ O.lcx (corresponging to 
SD ~ 0.005i?) over which the bridging transitions shown in Fig. [5] are smeared out. 

The thermodynamic states which are located on the metastable branches of the free 
energy curves survive during an average lifetime r rs tq exp^AQs/^BT) where AQ$ is the 
height of the energy barrier that separates the metastable from the stable branch and To is 
a characteristic microscopic time scale for the dynamics associated with the transition from 
a metastable to a stable wetting layer configuration. The energy barrier is highest in the 
vicinity of the bridging transition and vanishes near the ends of the metastable branches. 
An estimation of the energy barrier height yields, e.g., AQs ~ 75e for A/i = and D = 50a 
(a = 0.5R), and with ksT ~ e it follows that exp(AflsfkBT) ~ 10 32 , i.e., the metastable 
unbridged state for a = 0.5R near the bridging transition remains stable practically forever. 
However, at, e.g., a = 0.2R one has exp(Af2s//ceT) ~ 10 11 so that with r ~ lps. . . Ins one 
may observe a decay of the metastable states near the ends of the metastable branches within 
seconds or minutes. Thus the change of the morphology of the wetting films is expected to 
exhibit pronounced hysteresis effects as function of D. 

Obviously, in the limit of large separation D — > oo (in which only the configuration 
without a bridge is stable) the grand canonical potential Qs(-D) approaches the limiting 
value 20cP corresponding to the free energy of two individual spheres, each surrounded by 
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a wetting layer. It is convenient to separate this constant contribution 2tt s from the grand 
canonical potential n§ of the system and thus to define an excess free energy Qe(D) = 
Qg(D) — 2fig which contains all contributions from the wetting-layer induced interaction 
between the two spheres. In the limit D — > oo, i.e., in the absence of a liquid bridge, this 
excess free energy Qe(D) decays as D~ & (see Eq. ( |2.22| ) and, c.f., Sec. |V|). We note that 
for the example shown in Fig. |5| the coefficient of this leading order is positive, i.e., the 
effective potential in the absence of a liquid bridge is repulsive. This is owed to the choice 
T < T tt for this example: the spheres disfavor the adsorption of thick liquid films and the 
presence of the second sphere with its surrounding liquidlike layer leads to an additional cost 
in free energy which diminishes for increasing D. For the choice T > T u , i.e., if the spheres 
favor the adsorption of liquid (e.g., for T* = 1.3 as in Figs. |](a) and (b)) the coefficient 
of -D" 6 is negative and the effective interaction is attractive. However, in the presence of a 
liquid bridge, i.e., for sufficiently small values of D, the effective potential shows the same 
qualitative behavior as in Fig. |5| also for the case of thick wetting layers (T* = 1.3 > T* t ) as 
well as for the larger spheres (R = 50cr) with thin or thick films. 



B. Effective interaction potential for large spheres 

In this subsection we consider the limiting case that the sphere radius R is much larger 
than the diameter a of the solvent particles and that the separations a between the surfaces 
of the spheres are proportional to R: R^><7,a<^amR. For such large separations as com- 
pared to a the contributions Vt d (Eq. fl2.21|) ), Vt ei (Eq. (|2.16|) ), and $ (c.f., Eq. fl5J])) become 



vanishingly small relative to the contributions Qi g (Eq. (|2.18|) ) and Q ex (Eq. (|2.14|) ). For 
the case described above fli g and Q ex scale proportional to the surface area of the spheres, 
i.e., ~ R 2 , whereas for a/ a — > oo, R/a — > oo, a/R finite, $(-D; R) remains finite ~ e ss a & ss p 2 s 
with a proportionality constant of the order 1. Analogously, in the same limit Q e i — 2fl^ 
(Eq. ( [2.16D ) and Q s i — 20^] (Eq. ( |2.21| )) are determined by finite terms ~ Appiea 6 and 
~ App s e s fO-®f and of terms ~ pfea 6 and ~ pip s e s f(T®f, respectively, each with a proportion- 
ality constant of the order 1. Therefore measured in units of 87ri? 2 the unbridged branch of 
Qe in Fig. |5|(b) vanishes in the limit R — > oo. Moreover, on this scale the excluded volume 
at small a disappears from the figure, too, because d s /R — ► 0. 

Fig. P shows the excess effective interaction potential Qe in the limit of large spheres 
for the case A/i = 0, i.e., at two-phase coexistence in the solvent. In this limit and for 
A/i = 0, flig is the only relevant contribution to Qs because f2 ea; (A/x = 0) = 0. Accordingly, 
in this case the bridging transition is determined by the equality of the surface areas of the 
liquid-vapor interfaces for the unbridged and bridged configuration. From this condition 
and from dimensional analysis it follows that for large spheres D^(A// = 0) is determined 
by the equation 



8n(R + l Q fa\f = 8n(R + Z ) V£> / ( ^ ) (4.1) 



D 
R~+l 



where / is, for dimensional reasons, a universal function of D/(R+ Iq) alone which describes 
the surface area of the bridged configuration; Iq is the equilibrium wetting layer thickness 
on a single sphere. Since the line of bridging transitions lies below the line of thin-thick 
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transitions, Zo remains microscopicly small at the bridging transition (Fig. f|). Therefore one 
has 



£> w (A/i = 0) = AOR+Z ) (4-2) 

with a universal number 

A « 2.32 (4.3) 

determined by /(A) = 1 (compare Fig. If one applies this reasoning to Fig. one finds 
A ~ 2.39. Therefore even for R = 20a this macroscopic approximation leads to a surprisingly 
small error of only 3% for /^(A/i = 0). Accordingly, in Fig. |5] the full curves corresponding 
to A/i = closely resemble the ones in Fig. |6| describing the case of large spheres. The 
only differences appear for small separations a where the bridged branch linearly extends 
down to its minimum value Q E /{8n(R/a) 2 ) ~ — 0.0227e at a/R = (Fig. Only in this 
range of separations the effect of the contributions VL ei and Q s i becomes significant, leading 
to the deeper minimum visible in Fig. |5|. Thus for A/i = and large R the dependence 
of the effective interaction potential on R for the bridged configuration is captured by the 
indicated rescaling of the axes in Fig. |5](b) . However, our numerical analysis shows that the 
smallness of the deviations between the macroscopic description valid for R ^> a and the 
actual results for R = 20a is somewhat fortuitous. Whereas the dependence of Z)& f (A/x = 0) 
on R is indeed weak, the shape of the potential (for A/i = 0) reduces to that shown in Fig. |6| 
only for R larger than several hundred a and, surprisingly, for R up to 20 . . . 30<j with the 
deviations being maximal for R m 100a. 

Off coexistence AQb = &b(pi) — ^b(p g ) ~ A/tAp is positive so that Eq. (|4.1| ) has to be 
augmented correspondingly: 

8n(R + Z ) 2 + ((R + Zo) 3 - i? 3 ) = A + ^Vol(£) (4.4) 

3cr ig % 



where A and Vol(£) (Eq. Q2.15|) ) are the area of the liquid- vapor interface and the volume 
of the liquid, respectively, for the bridged configuration. They are obtained by inserting into 
Eqs. ( [2.19D and (|2.15|) that profile h(z) which solves the differential equation determining the 
minimum of Qi g [h] + Q ex [h] together with the appropriate boundary conditions. By splitting 
off a factor (R + Zo) 2 from A and (R + Z ) 3 from Vol(£) dimensional analysis shows that 
up to terms ~ Iq/R the critical distance for the bridging transition is given by a universal 
scaling function A: 

D U (A») = A ( ^] R (4.5) 



with A(0) = A. Thus off coexistence the critical bridging transition depends, apart from an 
explicit factor R, on R and A/i via the scaling variable ApA/iR/a^ . This property is shared 
by the whole bridged branch of the effective interaction potential. Thus increasing R for 
fixed undersaturation A/i has the same effect as increasing A/i for fixed R. From Fig. |^(b), 
in which the unbridged branch will disappear in the limit R^> a, one infers that the range 
and the depth of Qe decrease for increasing R at fixed undersaturation A/i. The behavior 
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of Dm and of the bridged branch of the effective interaction potential off coexistence and 
for R — > oo is determined by the behavior of the scaling function A(x) in the limit x — > oo. 
Our numerical data indicate that A(x — > oo) < 2 so that due to the geometric constraint 
D > 2R there is no bridging transition and the bridged branch of the effective potential 
vanishes for any value of A/x in the limit R —>■ oo. The cost in free energy due to the excess 
contribution Q ex suppresses the formation of a liquidlike bridge in the case of macroscopicly 
large spheres. In turn, this means that for any finite value of A/i there is a large but finite 
critical radius R c for which the critical separation a& t for the bridging transition attains the 
value a&t = 0, such that for R > R c there is no bridging transition. The determination of 
R c requires to analyze the full dependence of A on the scaling variable x. This, however, 
implies such a large numerical effort that it is beyond the scope of the present paper. 



V. DISCUSSION 

A. Total interaction potential 

The bare dispersion interaction between the two spheres is not included in Eq. ( [275] 
According to Hamaker |29| this contribution is given by 



, , A ss ( AR 2 AR 2 , (CD -2R)(D + 2R)\\ , , 

as the dispersion interaction between two identical spheres of radius R at center-of-mass 
distance D. In the limit a/R <C 1, where a = D — 2R (see Fig. |l](a)) is the smallest 
separation between the surfaces of the spheres, Eq. (|5.1|) reduces to 



A R 

§(D = 2R + a]R^>a)^—^-- (5.2) 

which corresponds to the Derjaguin approximation whereas Q(D ^> R; R) = 
— 16A SS R 6 /(9D 6 ). Thus except for the D-independent bulk contribution Vol(V 'f)^b(Pg) the 
total grand canonical potential of the system is 

{l tot (D; R) = {l s (D; R) + $(D; R) (5.3) 

where Qs{D]R) is given by the minimum value mxD.if l r x -\\(Tlg[h(z)]) for given D and R 
(Eqs. ( |2.13|) and ( |2.23| )) ; in analogy to Qe w e define the excess total free energy QE,tot = 
Qtot — 2Q^ . A ss is the Hamaker constant appertaining to the bare dispersion interaction 
between the particles in the spheres. In the case of pairwise additivity of the molecular 
interactions and in the absence of retardation effects one has A ss = An 2 € ss <j & ss p 2 if the inter- 
action potential between two individual molecules in the spheres is given by a Lennard- Jones 
potential (Eq. ( [2.4|) ) with the parameters e ss and a ss . Typically A ss is of the order of 10 _19 J 
or, equivalently, 10 . . . lOOe. If the vacuum between the spheres is replaced by a medium 
of condensed matter the interaction between the spheres is screened | 30| . In our present 



model this medium is the bulk vapor phase modified by the presence of the liquidlike films 
adsorbed on the spheres and the screening effect is described microscopicly by the functional 
Q[p(r)). 
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In Refs. |HJ and |32| this additional screening effect, due to spherical shells of adsorbed, 
homogeneous layers surrounding spherical particles, on the dispersion interaction between 
the latter immersed in another homogeneous medium has been calculated macroscopicly. 
Beyond molecular scales these results should closely correspond to the configuration without 
liquid bridge discussed herein because the deviation of the spherical shape of one wetting 
layer due to the presence of the second sphere is very small. Indeed, the interaction energy 
calculated in Refs. |TT| and is practically the same as the sum of the D-dependent 
contributions in fl ei (Eq. Q2.16Q ) and Q s i (Eq. ( |2. 21|) ) for the configurations without bridge 
- for these configurations Qi g and Q ex do not contribute to the dependence of Q$ on D - 
and the direct dispersion interaction $(-D; R). In Ref. [[H]] the total dispersion interaction is 
shown to be always attractive if the Hamaker constants corresponding to the interaction 
between any two media i and j are chosen such that Aij = An Ajj. Although the effective 
interaction induced by the wetting layers shown in Figs. ||](b) and |7| for the configuration 
without bridge is repulsive, we note that the sum Q tot of this interaction and of the bare 
dispersion potential $(-D; R) is also attractive if we choose the Hamaker constant in Eq. fl5.1|) 
accordingly, i.e., A ss = A 2 ^/A (Fig. [j]). Therefore our results are consistent with those 



obtained in Ref. ||31|| . Since only effective interactions between finite volumes enter into the 
total excess interaction potential flE,tot arid these effective interactions decay as D~ 6 in the 
limit of large separations D, the same holds for £lE,tot- 

Figure |] shows that as soon as the wetting films snap to a liquidlike bridge, whether it is 
stable or metastable, there is an attractive wetting-layer-induced force —dflE/dD that pulls 
the spheres together. From Fig. |5| one can infer that this attractive force is of the order of 
AOe/a in the range between a ~ 4<r (i.e., 0.2R for R = 20a discussed in this figure) and 
a ~ lOo" (0.5R) where the effective potential varies almost linearly. At the small separation 
a min ~ 2.5a the effective potential Qe induced by the bridgelike wetting layer is minimal and 
the wetting-induced force is zero. Finally, at smaller separations the interaction is repulsive 
leading to a stabilization of the spheres at D — D rnin = 2R + a m i n . Within the range a <C R 
the bare, direct dispersion interaction between the spheres (Eqs. (|5.1| ) and ( |5.2|) ) gives rise 
to a force F bare (a) « — A ss R/12a 2 . The estimate A ss ps 4ii 2 e 2 s j-aljp 2 s /ea 6 ~ 400e for the case 
of pairwise additive interactions without retardation follows from the ansatz A s f = \J A A ss , 
so that the bare dispersion force in our example with R = 20a is F& are (a) ~ — 670ea/a 2 . 
Therefore in the range where the bridge- induced force is almost constant (4cr < a < 10cr) 
the direct, bare dispersion force decays from approximately — 40e/a (which is of the same 
order of magnitude as the bridge-induced force) to approximately —6e/a, whereas for smaller 
separations it becomes the dominant force. 



B. Relevance for force microscopy 

Our model calculations can be tested experimentally by force microscopy. This can be 
done by suitably fixing one sphere in the fluid and by attaching the second one to the tip of 
a force microscope. Alternatively, both spheres can be positioned by optical tweezers and 
the force law can be inferred by monitoring optically their dynamics after switching off the 
tweezers. At separations between the spheres which are comparable with the diameter a 
of the solvent particles the actual effective interaction potential will exhibit an additional 



oscillatory contribution due to packing effects which decays exponentially on the scale a [33 



15 



In order to obtain these oscillations one would have to resort to density functional theories 
which are more sophisticated than the one in Eq. (|2.5| ). This, in turn, would make it 
much more difficult to obtain the bridgelike configuration, to map out the complete phase 
diagram, and to obtain results for large spheres. According to Subsec. [IV B| , for R^> a and 
at two-phase coexistence A/i = the bridging transition occurs at distances a which are 
proportional to R. In this case, due to R ^> a, the effective interaction potential will be 
practically unaffected by this oscillatory contribution for the vast portion a *C a <C a&t of 
the range of the effective interaction potential. 



C. Relevance for charge stabilized colloidal suspensions 



Whereas the kinds of experiments considered in the previous subsection are focused on 
two individual spherical particles, we discussed in the Introduction that the effective inter- 
action potential enters into the collective behavior of colloidal suspensions such that the 
bridging transition may trigger flocculation. If colloidal suspensions would be governed by 
dispersion forces alone, most of them would flocculate even in the absence of the wetting- 
induced forces discussed here because the dispersion forces generate the so-called primary 
minimum in the effective interaction potential close to contact. Since this minimum is much 
deeper than fcgT the colloidal particles would simply stick together permanently. This ef- 
fect, which is undesired for many applications, can be avoided by endowing the particles 
with electrical charges which adds a screened Coulomb repulsion between the charged par- 
ticles. As a result, such charge stabilized colloidal suspensions are characterized by effective 
interaction potentials in which a substantial energy barrier separates the aforementioned 
primary minimum from a second, much more shallow minimum at larger distances. Since 
this potential barrier is typically large compared with ksT the phase behavior of the col- 
loidal particles is practically independent of the primary minimum formed by the dispersion 
forces and determined by the shape of the potential outside the barrier. As demonstrated by 
Figs. [|and ^ the range of the wetting-induced forces is about 0.55.R, in good agreement with 
D bt ps 2.32(i? + l ) (see Eqs. fl4j) and fl4~3|)). On the other hand the position (and height) 
of the aforementioned energy barrier depends sensitively on the size of the total charge on 
the spheres, the amount of salt in the solvent, and the dispersion forces and can be varied 
over a wide range. With a high salt concentration the barrier position can be as small as a 
few tens of nm. Thus under such circumstances the wetting-induced interaction potentials 
would be relevant even for colloidal particles whose radii are only a few tens of nm. 



D. Relevance for stericly stabilized colloidal suspensions 



There is another class of colloidal suspensions for which the wetting-induced forces can be 
of practical importance. By coating the colloidal particles with polymers and by matching 
the refractive indices of the colloidal particles and the bulk fluid (in our case study the vapor 
phase or, more realisticly in the present context, the A-rich liquid phase of a binary liquid 
mixture acting as the solvent) the colloidal particles behave effectively like hard spheres 
(see, e.g., Refs. and [pSfl). Through this index matching the sum of the bare interaction 
potential §(D; R) and the effective interaction potential fl sg , which would arise if the spheres 
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were immersed in the homogeneous and unperturbed bulk solvent, vanishes. Within our 
model Q sg is given by the expression in Eq. ( j2.21| ) with p\ replaced by p g , which is the density 
of the bulk phase. Since the index matching works for the bulk phase, it does not work for the 
wetting phase. As a consequence the wetting-induced forces appear against a background 
effective potential of hard spheres. Therefore for this class of colloidal suspensions the 
wetting phenomena discussed here are expected to have a pronounced effect on their phase 
behavior. Within our model, for zndex-matched suspensions the total effective interaction 
potential is given by 

n totjjm {D- R) = n tot (D; R) - ($(£>; R) + il sg (D; R)) = Q S (D; R) - il sg (D; R) (5.4) 

and in analogy to Qe and ^E,tot we define 

n E ,im(D; R) = Q totiim (D; R) - 2fi^(i2) (5.5) 

with Q E>im (D — > oo; R) = for the unbridged solutions. Figure |8| displays VL tot ^ im and &E,im 
as function of a = D — 2R for the same system as in Figs. || and [7|. Q sg is about 30% smaller 
than Q s for the unbridged solution and also approaches its asymptote 2fiip from above. As 
before (see the discussion of Fig. |7] above), the resulting total effective interaction between 
spheres in an index-matched bulk fluid for the state with liquid bridge is still attractive and 
of the same order of magnitude as the bare dispersion interaction between the spheres, i.e., 
in the absence of the solvent. 



VI. SUMMARY 



We have obtained the following main results: 

1. Based on microscopic interaction potentials and within a simple version of density 
functional theory (Eqs. ( p.5|) -( ^9|) ) we have calculated the grand canonical potential 
of a system of two spheres immersed in a bulk fluid phase (Fig. p. The microscopic 
interactions are chosen such that the spheres prefer the adsorption of a second fluid 
phase which is thermodynamically close to the bulk fluid phase. Accordingly, a single 
sphere immersed in the fluid is covered by a homogeneous wetting layer of this second 
phase of thickness Iq. These thin wetting layers covering the spheres lead to an effective 
wetting-induced interaction potential fls(D) between the spheres. We have systemati- 
cally determined the dependence of Qs on the distance D between the spheres in terms 
of the morphology h(z) of the wetting film enclosing the spheres (Eqs. ( |2.13| )-( |2.21| )). 
We find that the shape of the effective interaction potential tts(D) depends, inter alia, 
on the effective interaction of two spheres immersed in the homogeneous wetting phase 
(Eq. ( |2.21| )). This contribution, which is independent of h(z), is not incorporated in 
previous phenomenological models for this system p2|.[2l 



2. The equilibrium interfacial profiles of the wetting layers are determined numerically 
by minimizing the free energy functional f2s[/i(z)] in Eqs. ( 2.13 )-( 2.21|) . We have cal- 



culated the rich structure of these equilibrium profiles (Fig. ^) for spheres of radii 
R = 20cx (Fig. H) and R = 50a (Fig. ||) where a denotes the diameter of the solvent 
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particles. As function of distance D, temperature T, and undersaturation A/i the 
system undergoes a first-order "bridging transition" between the two configurations 
shown in Fig. [T| For a fixed distance D we have mapped out the phase diagram of 
bridging transitions in the T-A/x plane (Fig. It turns out that the bridging transi- 
tion differs from and to a large extent is independent of the thin-thick transition of the 
wetting layer on each single sphere which is a remnant of the prewetting transition on 
the corresponding flat substrate. Thus one has to distinguish between the prewetting 
line for a first-order wetting transition on a planar substrate, the thin-thick transition 
line for wetting on a single sphere, and the bridging transition line for two spheres 
(Fig. ^|). At two-phase coexistence A/i = and for R ^> a the bridging transition is 
determined by the equality of the surface areas of the interfaces in the bridged and the 
unbridged configuration, leading to a universal ratio D&t(A/i = 0)/(R + lo) ~ 2.32 for 
the critical distance Db t (Afi = 0) of the bridging transition at coexistence (Fig. ^| and 
Subsec. |1V B| ). Off coexistence D bt (Afj,, R) is described by a universal scaling function 
(Eq. O). 



3. At large distances and depending on the temperature relative to the thin-thick transi- 
tion temperature on a single sphere the wetting-induced effective interaction potential 
can be either attractive or repulsive; in both cases it decays ~ D~ % for large D. The 
bridging transition leads to a strong break in slope of the effective interaction poten- 
tial at D = Dtt- This is the fingerprint of a cooperative phenomenon among the fluid 
particles whose degrees of freedom have been integrated out (see Sec. |). Metastable 
branches of the effective potential give rise to pronounced hysteresis effects (Fig. ^). 

4. In the case that a bridge of the wetting phase connects the spheres (i.e., D < Dm) there 
is an attractive wetting-induced interaction (Fig. ^) that pulls the spheres together. 
Within a wide range of separations a = D — 2R of the spherical surfaces this force is of 
the same order of magnitude as the bare dispersion interaction potential $ (Eq. (|5.1| )) 
between the spheres. This bare interaction of two spheres (corresponding to the case 
that they are separated by vacuum) has to be added to the effective potential Qs to 
yield the total interaction potential Q tot between the spheres which is attractive at 
large distances (Eqs. (|5.1| ), ( |5.3| ), and Fig. ^). 

5. The wetting- induced force between spherical particles is experimentally accessible di- 
rectly through suitable force microscopy techniques (Subsec. |V B| ). Moreover, in Sub- 
sec. |VD| we argue that this force influences the phase behavior of stericly stabilized, 



index-matched colloidal suspensions. The total effective interaction potential for such 
a case is shown in Fig. § it is repulsive at large distances. The phase behavior of charge 
stabilized colloidal suspensions (Subsec. |VC|) is only affected by the wetting-induced 



interaction potential if the screening length of the Coulomb repulsion in the solvent is 
smaller than a bt = D bt — 2R « 0.32_R. Depending on the size of the charges, the salt 
concentration of the solvent, and the underlying dispersion forces this criterion may 
be fulfilled even for colloidal particles whose radii are only a few tens of nm. 
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APPENDIX: CONTRIBUTIONS TO THE FREE ENERGY 



Our choice of interaction potentials <j>(f) (Eq. ( |2.1|) ) and <p s {r) (Eq. ( |2.4|) ) leads to the 
following expressions for the contributions to the free energy Qs (with the thermodynamic 
limit already carried out): 

oo oo 

n ei [h{z)] = 2Ap I Pl j dz{g + {z)+g_{z)) - J dz (f+(z) + /_(*)) | (Al) 





with 
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where R% = R + d s and 



4 V7V(^± + ^) 7 (k±-R) 7 J \(k± + R) 8 (k±-R) 8 

-™ 3 " k^R + R {jkTTW + JkT^w) ) (A3) 

with k± = ^h 2 (z) + (z ± D/2) 2 . The contribution is given by 

oo oo 

o \h( W 3/ A , f 1 / ^y° +yl(2g 2 +P 2 )+3p 2 gyi+p 4 \ 

w*)]=-w(A/o y^y^ — (^+2^+^)3/2 p ) 



1 ( qyj + y%{2q 2 + p 2 ) + 3p 2 qy 2 + + p 4 \\ 
%U (p 2 + 2y 2 +q + yXfl 2 P JJ [ } 

where the abbreviations y±, p, and q are defined by 

y 2 ± = a 2 + (z±z') 2 , (A5) 

p = h 2 (z) -h 2 (z'), (A6) 
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and 

q = h 2 {z)+h 2 {z'). (A7) 

The double integral in Eq. ( |A4| ) demonstrates the nonlocal functional dependence of Qi g on 
h(z). 
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FIGURES 





FIG. 1. Wetting film (thick full line) surrounding two identical homogeneous spheres of radius 
R which are separated by a distance D. The whole system is rotationally symmetric around the 
z axis which runs through both centers of mass. The position of the liquid- vapor interface which 
encloses both spheres is described by a function h(z), i.e., in cylindrical coordinates the sharp 
interface is given by the manifold {r(r±,(p, z) = (rj_ cos (p, rj_ sin (ft, z) G M 3 |rj_ = h(z)}. The origin 
of the coordinate system is in the middle between the two spheres so that their centers are located 
at z = ±D/2. a = D — 2R is the shortest separation between the surfaces of the spheres. Within 
the so-called sharp-kink approximation this interface separates a region of constant liquid number 
density pi from the surrounding bulk vapor phase of constant number density p g . Close to the 
surfaces of the spheres the repulsive interaction leads to a volume with thickness d s excluded for 
the centers of the fluid particles. For sufficiently large values of D the bridgelike wetting film 
configuration shown in (a) breaks up into two disjunct pieces so that h(z) = for a finite interval 
around z = (b). 
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FIG. 2. Morphologies of liquidlike wetting layers on two adjacent, identical spheres with radius 
R = 20a. The center-of-mass distance between them is D = 50a. The pictures show cross-sections 
through the system defined by the plane y = 0; the system is rotationally symmetric around 
the z axis (see Fig. |l]). The thick full lines denote the liquid- vapor interface, the thin dashed 
lines the surfaces of the spheres, (a) and (b): layer configuration with and without liquid bridge, 
respectively, for the temperature T* = ksT/e = 1.3 > T t * and at liquid-vapor coexistence Afj, = 0. 
Because of its higher free energy the configuration without bridge is metastable (c.f., Fig. Q). 
These configurations are characterized by the interaction potential parameters us = 6.283ecr 3 , 
ug = 0.838ea 9 , and d s = a. The temperature is above the thin-thick transition temperature 
T t * ~ 1.271 for each single sphere. In (c) the interaction parameters and D are the same, but 
the temperature T* = 1.2 is below the thin-thick transition temperature T^ t , so that the wetting 
layer around a single sphere is thinner than in (a) and (b). Also at this temperature the bridge 
configuration is the stable one (c.f., Fig. ||). 



24 




FIG. 3. Morphologies of liquidlike wetting layers on two adjacent, identical spheres with radius 
R = 50a, D = 120cr, and for the same choice of interaction parameters as in Fig. ||. The thick 
full lines denote the liquid-vapor interface, the thin dashed lines the surfaces of the spheres. These 
pictures magnify the region between the spheres. The temperature is T* = 1.2, which is above the 
thin-thick transition temperature T 4 * ~ 1.191 at coexistence for a single sphere, and the pictures 
differ with respect to the undersaturation: Afi* = A/i/e = 0.05 in (a), 0.015 in (b), and 0.01 in 
(c). Between (a) and (b), at A/ij^ ~ 0.0235 the system undergoes the first-order transition from 
the state without bridge to the state with bridge, and at A/z^ « 0.0103 between (b) and (c) there 
is a thin-thick transition of the wetting layer around the single spheres which is the remnant of the 
prewetting transition on the corresponding flat substrate (see, c.f., Fig. Q). Note that in (b) there 
are six turning points (•) of the profile h(z) whereas in (c) there are only two. 
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FIG. 4. Temperature-undersaturation phase diagram of wetting layer configurations for two 
spheres with R = 20a at a fixed distance D = 50a (a = 10a). The interaction potential parameters 
are the same as in Fig. |2[ The three configurations shown in Fig. ||(a)-(c) are located at the 
respective thermodynamic states "a" to "c" (♦). The line of liquid- vapor coexistence A/i = 
separates the region where in the bulk the vapor phase is stable and the liquid phase is metastable 
(—A/i < 0) from the region where the liquid phase is stable and the vapor phase is metastable 
(—A/i > 0). The dotted "met ast ability line" ("ml") separates the region where the liquid phase 
in the bulk is still metastable (—A// > — Afj, m i(T)) from the region where only the vapor phase is 
stable in the bulk (—A/i < — A/i m ;(T)). The liquidlike layer on each individual sphere exhibits a 
first-order thin-thick transition at —A/i = — A/^(T) (dashed line "tt"). This line intersects the 
liquid-vapor coexistence line at T t * ~ 1.271 and ends at a critical point (•) in the vapor phase 
region; T t * c ~ 1.275 and — A/Xj i c ~ —0.0144. For the present choice of interaction potential 
parameters, at lower temperatures and larger undersaturations —A/i = — A/i^T) (full line "bt") 
the first-order bridging transitions between the configurations with bridge (— A/i > — A/i^(T)) 
and without bridge (—A/i < — A/i^(T)) occur. This line intersects the coexistence line linearly. 
Within the sharp-kink approximation the line of bridging transitions happens to be cut off by the 
"met ast ability line"; within a more sophisticated approach the line "bt" is expected to end at a 
critical point, too. The locations of the thin-thick transitions in the phase diagram are practically 
not affected by the presence of the bridge. The dashed-double-dotted lines (—••—) are metastable 
extensions of the thin-thick and bridging transition lines, respectively. The dashed-dotted line "p" 
(— • — ) is the prewetting line for the corresponding planar substrate. It joins the liquid-vapor 
coexistence line A/i = tangentially at the first-order wetting transition temperature ~ 1.053 
(A) and ends at a critical point (■) in the vapor phase region. For a discussion of the effects of 
fluctuations on this mean- field phase diagram see the main text. 
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FIG. 5. (a) Dependence of the grand canonical potential fig on the separation a = D — 2R 
and the undersaturation Afi for the same system as in Figs. ||(c) and ^ i.e., for R = 20a and 
T* = 1.2 < TZ. The dots indicate the end points of metastable branches. For a < 2d s (with 
d s = a = 0.05i? here) the excluded volumes around the spheres overlap. In the limit D —* oo 
the stable solution is the one without a liquid bridge; in this limit fis(-D —* oo) = 2D,g is twice 
the free energy of a single sphere surrounded by a wetting layer. At the separation D^t or c%, 
where the two free energy branches intersect for a given Afj,, a first-order morphological phase 
transition between a configuration with a liquid bridge and a state without bridge takes place. The 
equilibrium thickness of the homogeneous wetting layer around a single sphere is Iq ~ 1.3a, so that 
Df,t/(R + Iq) ~ 2.39; the slight deviation from the prediction of Eqs. (O) and ( p~3| ) is due to the 
still rather small size of the spheres. We note that, in contrast to the case shown here, for T > Tu 
the free energy curve corresponding to the solutions without bridge approaches its asymptote from 
below, (b) Same as in (a), showing the excess free energy £Ie = &>s ~ 20& . In this presentation 
the results for the solutions without bridge and for different undersaturations A/x collapse onto a 
single line. Qe{D — > oo) decays as D~ 6 . 
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FIG. 6. Excess free energy Ue = ~ 2fi<j for A/i = in the limit of large spheres, i.e., 
R 3> cr, a -C a ~ R. In this limit the excess free energy branch for the unbridged solution vanishes 
if it is measured in units of 8nR 2 . Off two-phase coexistence, i.e., for Ap ^ the branch for the 
bridged solution is determined only by the contributions (Eq. ( |2.18| )) and £l ex (Eq. ( |2.14 )) to 
the free energy. At two-phase coexistence A// and Sl ex vanish so that Qe is solely determined by 
Qi g . Therefore within the local theory with ^f"^ (Eq. Q2.19| )) the bridged solution is a minimal 
area surface, i.e., its mean curvature is zero. Since a 3> d s ~ a the excluded volume at small 
a disappears from the figure. Therefore, compared with the full curve in Fig. |5|(b) the potential 
curve here is effectively shifted to smaller values of a. Moreover, the actual minimum of the 
effective interaction potential at small a > o (compare Fig. ||), which is due to the influence of 
the contributions f2 e j and £l s i, is not visible on this scale either. The critical separation for the 
bridging transition (♦) is given by a^t/R ~ 0.32 (Eqs. (O) and (|4.3j )). If the thermodynamic state 
of the system is driven into the off-coexistence region A/i > the whole excess free energy branch 
for the bridged solution is shifted upwards (compare Fig. |B|). For any finite value of A/x, in the 
limit R — » oo there is no bridging transition anymore (see the main text). 
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FIG. 7. Excess free energy Qe = ~ 2Q S (dashed lines) and excess total free energy 
Q E tot = n s - + $ (full lines) for A/i = 0. Here T* = 1.2 and R = 20a so that the 

dashed lines are identical with the full lines in Fig. ||(b). The dots indicate the end points of 
metastable branches. The parameters e ss and a ss of the pair potential between the particles 
forming the spheres are chosen such that the condition A s f = ^ AA SS for the corresponding 
Hamaker constants is satisfied. Although the wetting-layer induced potential for the solutions 
without bridge is repulsive, the total interaction potential including the bare dispersion potential 
is attractive. For small separations a or D the bare dispersion potential dominates. In the limit 
D — > oo, i.e., for the configurations without bridge, Qe and &>E,tot decay as D~ e as expected for 
dispersion interactions. 
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FIG. 8. Same as in Fig. || but with Q t ot,im = ~ ^sg (a) and with 0,E,im = ^tot,im — 2 ^!m 
(b). We choose again T* = 1.2, R = 20a, and the interaction parameters as in the previous 
figures. The dots indicate the end points of metastable branches. The total interaction potential 
for index- matched spheres and bulk fluid is again repulsive: since the temperature is below the 
thin-thick transition temperature Tu the adjacent spheres dislike the presence of additional liquid in 
their vicinity and therefore it is energetically advantageous to separate them as much as possible. 
^E,im for the solutions without bridge is smaller than f2#. However, for the bridged solutions, 
&E,im and £Ie as well as the corresponding wetting-induced forces are of almost the same size, 
respectively. 
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